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ABSTRACT 

We compute the B-mode polarization power spectrum of the CMB from an epoch of inhomogeneous 
reionization, using a simple model in which H II regions are represented by ionized spherical bubbles 
with a log normal distribution of sizes whose clustering properties are determined by large-scale 
structure. Both the global ionization fraction and the characteristic radius of H II regions are allowed 
to be free functions of redshift. Models that would produce substantial contamination to degree scale 
gravitational wave B-mode measurements have power that is dominated by the shot noise of the 
bubbles. Rare bubbles of > 100 Mpc at z > 20 can produce signals that in fact exceed the B-modes 
from gravitational lensing and are comparable to the maximal allowed signal of gravitational waves 
(~ 0.1 /jK) while still being consistent with global constraints on the total optical depth. Even bubbles 
down to 20 Mpc at z ~ 15, or 40 Mpc at z ~ 10 can be relevant (0.01 (iK) once the lensing signal 
is removed either statistically or directly. However, currently favored theoretical models that have 
ionization bubbles that only grow to such sizes at the very end of a fairly prompt and late reionization 
produce signals which are at most at these levels. 

Subject headings: cosmic microwave background — cosmology: theory — large-scale structure of 
universe 



1. INTRODUCTION 

The polarization of the cosmic microwave background 
(CMB) can be decomposed into two components, E- 
modes and B-modes, which are geometrically distinct 
but non-local combinations of the Stokes parameters Q 
and U . The two modes are distinguished by the rela- 
tionship between the direction of polarization and its 
spatial gradients. In linear perturbation theory, both 
scalar and tensor perturbations contribute to the E- 
modes, but scalar perturbations do not produce B- 
mode s (|Kamionkowski et alj|1997t iSeliak fe Zaldarriagal 
Il997f h This fact makes observation of the CMB J5-modes 
an ideal way to detect the tensor perturbations from 
gravitational waves generated during inflation. Detection 
of gravitational waves would be strong evidence in favor 
of inflation, would determine the energy scale of inflation, 
and would help constrain models of infl ation (see e.g. 
iGold fc Albrechll200a IVerde et alJl200l for recent as- 
sessments). This detection is the goal of a concerted ex - 
perimental effort currently underway (jBock et alJ [2006). 

Beyond linear perturbations, scalar density fluctua- 
tions can generate B-modes. For example, gravitational 
lensing spatially distorts the polarization field in a man- 
ner that is independent of its orientation and hence gen- 
erates B- modes from the intrinsic B- modes from recom- 
bination (Zal darriaga fc Sel iak 1998). Lensing places a 
fundamental limitation on the detection of inflationary 
gravitational waves that is not far in energy scale from 
current limits (e.g.. [Page et al.l l2006f ) if only polariza- 
tion power spectra are measured. Fortunately, with high 
signal-to-noise, high resolution polarization maps, de- 
lensing of the polarization field is in princip le possible 
(|Hu fc Okamotol 120021 : iHirata fc SeliaM I2003D and may 
improve limits by an order of magnitude or more in po- 



larization power CKnox fc SonJ2002HKesden et all20(r2i 
ISeliak fc Hiratal l2004h . In addition, predictions for the 
gravitational lensing B-modes are highly accurate and 
hence their contamination to a power spectrum measure- 
ment may be subtracted to within a few percent for a 
cosmic variance limited experiment. 

More worrying would be cosmological B-mode sources 
that cannot be accuratel y modelled theoretically. As 
pointed out in iHul (|2000| ). inhomogeneities in the free 
electron density during reionization will generate B- 
modes through Thomson scattering. Fluctuations in the 
local ionization fraction during an extended epoch of 
reionization are in fact expected. Current observations 
show evidence for an extended epoch of reioniz ation that 
ends at redshift z ~ 6 and begins at z > 10 (Fa n et aTl 
2006). Because reionization is a gradual process and 
the sources of ionizing radiation are clustered, the dis- 
tribution of neutral and ionized hydrogen during the 
epoch of reio nization is expected to be highly inhomo- 
geneous fe.g- lBarkana fc Loeb|[200ll ). 

Previous studies of B-modes generated by inhomoge- 
neous (or "patchy") reionization have shown that they 
are negligible compared with B-modes from lensing at 
scales relevant for gravitational waves if the ionized re- 
gions are small and un correlated with each o ther and 
reionization is prompt (|Hul 2000; Wcller fl999f ). Corre- 
lations between io nized regions can boost the signal but 
not substantially (iKnox et al.l 119981 : ISantos et all 120031 : 
iBaumann et al.lf2003f ). However, because reionization is 
a highly uncertain process, it has not been possible to 
rule it out as a source of significant B-modes when these 
conditions are violated. 

In this paper, we take a phenomenological approach 
to the study of B-mode generation from inhomogeneous 
reionization to bound the possible contamination to grav- 
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itational wave studies. We limit our considerations to 
models that satisfy observatio nal bounds on the total 
optical depth r» = 0.09 ± 0.03 (jSpergel et al.ll2006h . We 
fix the other cosmological parameters to also be consis- 
tent with the 3-year WMAP results: Hubble constant 
Hq = lOO/i km s _1 Mpc -1 with h = 0.73, matter density 
Q m h 2 = 0.f28, baryon density Q b h 2 = 0.0223, flat spa- 
tial geometry = 0, CMB temperature Tcmb = 2.725 
K, primordial helium fraction Y p — 0.24, scalar index 
ns = 0.958, and amplitude <jg — 0.76. 

Our approach should be contrasted with studies of 
physically motivated m odels of reionization (jLiu et alJ 
12001b IZahn et ai1l2005l ). The advantage of our approach 
is its ability to describe a wide variety of reionization sce- 
narios with relatively few phenomenological parameters. 
The parameters that are required to generate substantial 
-3-modes can then be compared with predictions from 
physical models. 

The outline of the paper is as follows. In Section [21 we 
introduce a phenomenological model of inhomogeneous 
reionization that incorporates several free functions of 
rcdshift. We parameterize these functions in Section [3] in 
a manner which highlights the phenomenology of the B- 
modes discussed in Section H) We discuss observational 
and theoretical bounds on the parameter space in Sec- 
tion [5] and conclude in Section [6l 

2. INHOMOGENEOUS REIONIZATION MODELS 

Reionization generates polarization by Thomson scat- 
tering of anisotropic radiation coming from recombina- 
tion. If the ionized hydrogen is homogeneous, the scalar 
quadrupole temperature anisotropy creates only .E-mode 
polarization. Inhomogcncities break this symmetry and 
generate -B-mode polarization. Under the Limber ap- 
proximation, valid for multipoles £ 3> 10, the 5-mode 
power spectrum is related to th e power s pectrum of the 
ionized hydrogen Ps HII (k, z) as (|Hul l2000) 
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where <jt is the Thomson scattering cross section, n p ^ 
is the present number density of protons, Zi and Zf are 
the redshifts at the beginning and end of the epoch of 
reionization, H(z) is the Hubble parameter, Da(z) is 
the angular diameter distance, r(z) is the optical depth 
out to redshift z, and Q 2 ms (z) is the variance of the 
quadrupole. We define the ionization fraction to equal 
1 when both hydrogen and helium are fully ionized, so 
n p> o — (l~^ / p/2)nb,o- In the fiducial model Q mls ~ 18 fiK 
at the redshifts of interest. We often use the notation 
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for the power per logarithmic interval in multipolc. 

To evaluate equation (TTJ) , we need a model for inhomo- 
geneous reionization. Because we are interested in ex- 
ploring what types of ionization histories generate max- 
imal amounts of -B-mode s, we adopt the phe nomenolog- 
ical parameterization of IWang &; Hul (120061) instead of 
a mo re physically motivated one fe.g.. IFurlanetto et al.l 
l200l . 

In this parameterization, each ionized H II region is a 
spherical bubble of radius R with full ionization inside, 



and neutral h ydrogen outs i de. W e generalize the param- 
eterization of I Wang &: Hul (|2006l ) to describe the power 
spectrum of H II regions rather than H I regions and al- 
low the H II bubbles to have a distribution of radii, P(R), 
that can evolve with redshift. Note that we suppress the 
redshift arguments to functions such as P(R) where no 
confusion will arise. The physical mechanisms of reion- 
ization, including what the ionizing sources are and how 
efficiently they can ionize neutral gas surrounding them, 
are left unspecified. Therefore not all of the parameter 
space available corresponds to currently favored, or even 
physically plausible, models. 

In Section 12.11 we describe the basic features of the 
model and the relationship between the bubble distribu- 
tion and mean ionization. In Section 12.21 we show how 
the ionization history and bubble radius distribution de- 
termines the power spectrum of the ionized hydrogen. 
Finally in Section [2.3l we explore the impact of the width 
of the distribution on the power spectrum. 

2.1. Mean Ionization 

We assume that H II bubbles populate a large 
volume Vq as a Po isson process with fluct uating 
mean (IValageas et al.l 120011: IFurlanetto et alJ 12004 
IWang fc Hd 12006ft . In addition we allow the bubbles to 
have an arbitrary distribution of radii P{R)- The volume 
of each bubble is Vb(R) = AttR 3 /3. Then the ionization 
fraction is 



(*(r))i 



1 — exp 



JdR^V b (R) 



(3) 



where dn^/dR — nb(v)P(R) and rib(r) is the comoving 
bubble number density, and the subscript P denotes av- 
eraging over the Poisson process. 
The average ionization fraction in the volume Vq is 

1 



1 



d :i 



v»j Vo -^[~J dR ii Vb{R) Y (4) 

We assume the number density of bubbles fluctuates as 
a biased tracer of the large scale structure with bubble 
bias b{R): 



drib _ dm r-. 
~dR ~ ~dR [ 



b{R)6 w (r)}, 



(5) 



where 8w{y) = J d 3 r' ' 5(y')Wr(y — r') is the overdensity 
6 smoothed with a top-hat window function 

W R (k) = ^3 [sin(fci?) - kRcos(kR)}. (6) 

In our model, we assume that the bubble bias is inde- 
pendent of R, b(R) = b for simplicity. 

Using equation ([5]) and expanding the exponential in 
equation (j4|), we find that the average number density, 
bubble volume and ionization fraction are related by 

p — n b (V b ) f f f Af. 

! I dR^V b (R)S w (r) 



1 



V 



d 3 r 



V'o 



= e -n b {V b ) 



where we define 



dR P(R)V b (R), 



(7) 



(8) 



and assume J v d^rSwi?) = for a sufficiently large vol- 
ume. Therefore any two of {x e , fib, (Vb)} specifies the 
mean ionization model. 
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2.2. Ionized Hydrogen Power Spectrum 

The ionization field can be written as a;(r) = x e +Sx(r), 
where x e is the spatially averaged ionization fraction of 
equation ([7]). The ionized hydrogen perturbations have 
contributions from both ionization fluctuations 5x and 
hydrogen overdensities 5 ()Hul [20001 

<W r ) = Sx(t) + [x e + 8x(t)]S(t). (9) 

The 2-point correlation function of ionized hydrogen, 
&w<w(r) = (^Hii(ri) JHii(r 2 )) where r = |r x - r 2 |, is 
(|Furlanetto et al.ll2004h 

& H nfen ( r ) = &x8x (r) + x 2 e £ss (r) (10) 
+£,Sx8 5x5 (r) + 2x e £, Sx $(r), 
neglecting terms that are higher order in S and Sx (except 
Csxd 6x5 1 which turns out to be important to this order 
in the perturbations). The H II power spectrum is the 
Fourier transform of the correlation function, 



P6mx(k) = J d 



' 3 r e ik ' 



^HII^HII v")' 



(11) 



Following Wan g fc Hul (|2006f ). we compute the correla- 
tion functions in equation (|10|) by separating the con- 
tributions from one-bubble and t wo-bubble correlations , 
analogous to the halo model (see lCoorav fc Shethll2002l 
for a review). 

2.2.1. Two-bubble terms 

For the two-bubble contribution to the two-point 
H II correlation function, the term £|£ a Sx5 is 
higher o rder than the other terms and can be 
dropped (jFurlanetto et al J 12004). We choose to include 
the density correlation term x 2 e ^ss with the two-bubble 
terms due to their joint dependence on the matter power 
spectrum P(k). The remaining two-bubble terms in the 
H II power spectrum are Pi xSx {k) and P 2 xS (k). To find 
these terms, we need an expression for Sx(r) that ac- 
counts for the clustering of H II regions through the bub- 
ble bias b. 

Expanding the exponential in the Poisson-averaged 
ionization fraction equation §3§ gives 



1-6 dR n b P(R)V b {R)5 w (r) 



(x{r)} P =l- e -"oW 



= x e + (1 - x e )fi b b J dR P(R)V b (R)5 w (r), (12) 

so the ionization fraction perturbation is 

(fa(r)) P = (z(r)) P - x e (13) 
(1 - a; e )ln(l - x e ) 



(Vh) 



b / dR P(R)V b (R)S w (r). 



Using this expression for Sx, we can find the correla- 
tions of Sx with Sx, and Sx with S. The Fourier trans- 
forms of these correlation functions are 

Plx5x(k) = [(1 - x e ) ln(l - x e )b(W R (k))} 2 P{k), (14) 

P? xS {k) = -(1 - x e ) ln(l - x e )b(W R (k))P(k), (15) 

where P(k) = Pss(k) is the matter power spectrum, and 
the two-bubble window function averaged over the bub- 
ble radius distribution is defined as 

(W R (k)) = [dR V b {R)P{R)W R {k), (16) 

\Vb) J 



where W R (k) is the Fourier transform of the real-space 
top- hat window of equation |6|). 

The total two-bubble contribution to the power spec- 
trum comes from combining the terms in equations (| 14[) 
and (fT5|) with the density fluctuation term P(k), with 
coefficients as given by equation (|10| : 



PsL( k ) = [(1 - *e)ln(l - x e )b{W R {k)) - x e ] z P{k). 

(17) 

Finally to complete the description we need a model 
for the bubble bias b. We assume that the bubbles are 
centered on dark matter halos so that the bias of the 
central halo determines the bubble bias. Specifically, we 
match the desired number density of H II bubbles to the 
number density of halos above a mass threshold M t h 



n b 



dM drih 



(18) 



where dnh jd In M is the halo mass function. In a physical 
model M t h might be associated with the threshold mass 
for collapse determined by coo l ing v ia atomic hydrogen. 
We use the iSheth fc Torrnenl (pflflh mass function for 
dnh / din M. The bubbles are then correlated according 
to the matter power spectrum but with a bias, 



h 1 r dM h <m\ dnh 

6 = t- / -rr b h( M ): 



n b Jm, 



M 



d\nM' 



(19) 



where the bias o f dark matter halos bh is given by the 
ISheth fc Torrnenl (|1999f ) form. Note that this association 
with dark matter halos is only necessary for modelling 
the two bubble correlations. For much of the parameter 
space of interest, the one halo term dominates and its 
contributions are independent of this association. 

2.2.2. One-bubble terms 

In the one-bubble regime, we can neglect in equa- 
tion (|10[) because we assume complete ioni zation inside 
bubb les and complete neutrality outside (|Wang fc Hul 
2006). The £55 term in £<s HII( s HII is already included in 
the two-bubble contribution to the power spectrum, so 
only the first and third terms in equation (|10[) remain. 
For our purposes, the main con tribution from gxS 
can be approxima ted by ^Sx^ss (jFurlanetto et al.ll2004l : 
I Wang fc Hij[200l . 

The two-point correlation of the total ionization 
fraction in the one- bubble regime can be written 
(jGruzinov 



(x(r 1 )x(r 2 ))=x 2 e +esx5x(r) 

= x e P 2ll {r)+x 2 e (l 



(20) 



P2\i(r)), 



where f2|i( r ) is th e probability of the point r 2 being 
ionized (in the same bubble as ri) given that ri is ion- 
ized, and r = |ri — r 2 |. The first term is the probability 
that both points are ionized in the same bubble, and the 
second is the probability that the points are ionized in 
separate bubbles. Then 



a|i 



dRP{R)V b {R)f (-^) , (21) 



where f(r/R) is the probability of ionization at r 2 given 
ri in an ionization bubble of radius R, with r — \y\ — r 2 |. 
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Suppose the point ri is in an ionized bubble of radius 
R at a distance ro < R from the center, and Y2 is sep- 
arated from ri by r. If we consider two spheres with 
radii R and r, with their centers a distance ro apart, 
then we can define A(R, r, ro) as the area on the sphere 
with radius r that is in the interior of the sphere with 
radius R. The probability that a random point ri inside 
a bubble of radius R is a distance ro from the center is 
P(r )dr Q = (ATtrldr )/V b {R). The function f(r/R) is 
the integral over ro from to R of the fractional area 
A(R, r, ro)/(47rr 2 ) times P(ro), which gives 

and f(r/R) — for r > 2_R. Assuming there is not 
significant overlap between bubbles, this function f(r/R) 
is the convolution of V b W R {r) with itself and its Fourier 
transform is V b W^{k) (|Wang fc Hull2006l ): 



f( 



\R 



l — J d\'V 2 (R)W R (v-v')W R {v') 



V b {R) 
V b (R) 



d A k 



W 2 R {k)e 



— zk-j 



(23) 



The overlap between bubbles will be small if n b (V b ) <C 1, 
which corresponds to x e -C 0.63 by equation ([7]). This 
assumption should be valid until near the end of reion- 
ization. 

Using this expression for f(r/R) with equations (|20[) 
and (1211). the ionization correlation function is 



£sx8x( r ) — [ X e 



<)P2\i(r), 



(24) 



; (fc)e 



— zk-: 



p ^ r) = W)S dRP ^ 2{R) Swf w ' 

We can take the Fourier transforms of CgxSx an d QxSx^ss 
to get the 1-bubble terms in the power spectrum 



8x8: 



(v b ) 



dR V 2 (R)P(R)W 2 R (k), (25) 



P} b x55x5 {k) = {x e -xi)P{k). (26) 

The first term, equation (|25|) . is associated with the 
shot noise of the bubbles. In the second term, equa- 
tion (f2"6"|) . P(k) is the convolution of the matter power 
spectrum P(k) with the square of the one-bubble aver- 
aged window function, 

P{k) = {V b ) J ^(wK^Pdk-kil), (27) 

where the one-bubble window function averaged over the 
bubble radius distribution is defined by 

(W R (k)) = j^JdR V 2 (R)P(R)W 2 R (k). (28) 



As in lWang fc Hul (|2006l ). we simplify our calculations by 
approximating P(k) by interpolating between its small 
and large k limits: 

P(k)(V b )(a R ) 



P{k) 



(29) 



[(P(fc)) 2 + (<H)K)) 2 ]V2' 

where (o~ 2 R ) is the top-hat smoothed density variance av- 
eraged over the bubble radius distribution 

(4) = ^JdR V 2 (R)P(R)a 2 R . (30) 



The total one-bubble contribution to the H II power 
spectrum is the sum of equations (|25|) and (|26|) . 



PsLW = a^ 1 - 1.) (H)<W3(fc)) + P(fc) 



2.3. Bubble Radius Distribution 
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To evaluate the power spectrum, we need to specify 
the bubble radius distribution which determines the av- 
eraged window functions (W R ( k)) and (Wpjk)). Moti- 
vated by the analytic mod el of iFurlanetto et al.l (|2004l ) 
and numerical simulations (|Zahn et al.l 2007f ). we assume 
here a log normal distribution with width <Ti n _R: 



P(R) = 



,-[ln(R/R)] 2 /2al 



(32) 



' In R 



This distribution is shown in Figure[T]for a\ nR — 0.2 and 
o-\ n R = 1.0. 

The one and two bubble terms in the power spectrum 
possess different weights over this distribution. From 
equation (|16p using a log normal P(R), the two-bubble 
averaged window function is 

(W R (k)) = ^= / d XX 2 e~ hi x/2 ^» (33) 



sin nx ~ K X cos K X 

(«x) 3 

where k = kR and x = R/R- Similarly, for the one- 
bubble terms, equation (|28|) becomes 



(W 2 R (k))- 



9e 



-9(T,, 



27TCTl n _R. JO 

(sin kx — kx cos «x) 2 



d X X 5 e~ ln2x/2<R (34) 



For a given value of a\ n R , each of these averaged window 
functions depends only on k = kR and not on k or R 
separately. Figures and [3] show the window functions 
for various choices of o~\ nR , evaluated by numerically in- 
tegrating the expressions in equations ([33]) and (|34|) . 

It is useful to examine the limiting behavior of these 
window functions to determine the characteristic scale 
in the average. Using the fact that lim^^o W R (k) = 1, 
we see from equation (|16p that the two-bubble window 
function (Wr) also approaches unity in the limit of small 
fc. In the oscillatory high fc regime, the integration will 
sharply suppress contributions. Thus the weight in the 
distribution simply reflects the average bubble volume 
and it is useful to define a volume weighted radius R 



which is larger than R (see Figure [T]) 
Ro = e 3 <n /2 R. 



(35) 



(36) 

The fc — * limit of equation (|28|) for the one-bubble 
window function is 



hm(^(fc)) 



(Vb) 



e 9 <«, 



(37) 



where the second equality holds for the case where P{R) 
is the log normal distribution. For large fc, W R (k) w 
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Fig. 1. — H II bubble radius distributions for (T\ n R = 0.2 
{dashed) and a\ n R = 1.0 (solid) with the same effective radius, 
Reft = 55 Mpc. The vertical lines show the location of the volume 
weighted radius Ro for each distribution. 




Fig. 2. — Squared two-bubble averaged window functions for log 
normal bubble radius distributions with widths cri n ji = 0.2 (thin), 
0.5 (medium), and 1.0 (thick). The square of a top-hat window 
function is plotted for comparison (dashed). 



9(kR) 4 cos 2 (kR). Since cos 2 (kR) oscillates rapidly in 
this limit compared with the variation of V 2 P(R)R~ 4 , 
we can replace (sin k\—kx cos k X) 2 by 1/2 in the integral 
in equation (|34p . Then we have 



lim (W 2 R (k)) : 



9e 



-9o,, 



■ In 2 X /2of, 



(38) 



To estimate the characteristic scale in the one-bubble 
terms, we take the value of k where (W R ) transitions 
between these two limits. This wavenumber, k c g, can be 




Fig. 3. — Same as Figure [2] but for one-bubble averaged window 
functions. 



estimated by requiring that 



}im(W 2 R (k)) 



lim (W 2 R (k)) 



Using equations ([37]) and ([38]) . we find 

1/4 



fc=fc r 



«eff 



R^e~ 2b < 



(39) 



(40) 



Therefore, we expect the main contributions to arise from 
scales 



RcB = Re 4 



Roe 



2..W 



(41) 



Note that the two distributions in Figure[T]have the same 
i? e ff even though their R differ by a factor of ~ 50. In 
the (Tin r —* limit, i? c ff = R = Ro so all the polarization 
power comes from bubbles with radii equal to i? c ff- For 
larger values of a\ n R, we find that most of the power is 
still due to bubbles with R ~ i? cff . 

3. PARAMETERIZATION OF REIONIZATION MODELS 

The framework for inhomogeneous reionization de- 
scribed in the previous section requires the specifica- 
tion of three functions of redshift: the ionization fraction 
x e (z), the volume weighted bubble radius Ro(z) defined 
by equation (|35() , and the width of the log normal radius 
distribution a\ n R(z). In this section we discuss several 
parameterizations of these functions that highlight their 
impact on the -B-modes. 

We begin by discussing the leading order shot noise 
term contributing to the 5-mode spectrum and how it 
determines the shape of the spectrum in Section |3~T1 We 
parameterize the range of redshifts for which < x e < 1 
by Zi and zt, the redshift at the beginning and end of 
reionization respectively. In Section 13.21 we explore the 
impact of various choices for x e (z) during this epoch, 
including histories that maximize the i?-modes and inte- 
gral constraints from the total optical depth . Finally 
we take two representative forms for Ro{z) in Section ^. 31 
a constant radius set by R c g the characteristic radius of 
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Fig. 4. — An example of a B-mode polarization spectrum from 
inhomogeneous reionization (solid, thick), and the contributions 
to the spectr um from each component of the H II power spectrum 
in equation 142 1 (thin, solid curves as labelled): the one-bubble 
shot noise term, the one-bubble convolution term, and the two- 
bubble term. The reionization model that produces this polariza- 
tion (model 19 in Table [2} has a linear ionization history, constant 
n b , z f = 6, Zi = 18.4, r* = 0.12, a lnR = 1.0, and R (z f ) = 50 Mpc 
(R c ff(zf) = 609 Mpc). The lensing B-modes for the assumed cos- 
mology (dashed) and the gravitational wave B-modes with tensor- 
to-scalar ratio r = 0.5 (dotted) are plotted for comparison. 



the bubbles for the dominant shot noise term, and a con- 
stant number density of bubbles which grow in radius. 
For simplicity we take a constant a^R throughout. 

3.1. Approximate B-mode Spectrum 

Different components of the H II power spectrum can 
vary greatly in how much they contribute to Cf at a 
particular scale. We study the relative contributions by 
separating the power spectrum into the one-bubble shot 
noise term, the one-bubble convolution term, and the 
two-bubble term 



P 



8mj 



plb , plfa 
r 8x8x ' r ! 



8x8 8x8 



p26 

^HII 



(42) 



and integrating equation {T]) replacing Ps HI1 by each of 
these three terms. Figure U shows these contributions 
to the total B-mode spectrum for a particular reioniza- 
tion model. The lensing and gravitational wave B-modcs 
plotted in Figure [5] show that B-modes from inhomoge- 
neous reionization would be an important contaminant 
for gravitational waves if they have a large amplitude at 
£ < 100, so we sometimes focus on the B-mode amplitude 
at £ ~ 100. 

For all values of model parameters we consider, the 
contribution to Cf from the convolution term, P} xS SxS , 
is less than the contributions from the other two terms. 

For some choices of parameters, the contribution from 
P} h is dominant over a range of £ values. In these cases, 
however, the B-mode amplitude is never larger than a 
few tens of nK 2 . The amplitude of the two-bubble con- 
tribution is limited by the fact that (Wi(k)) is never 
greater than 1, and because P(k) decreases with redshift. 



All of the models with large B- modes (Cf > 100 nK 2 ) 
at £ ~ 100 are dominated by the one-bubble shot noise 
term in the H II power spectrum, P^ Sx (e.g., the model 
in Figure [4]) . Note that this is the only term in the H II 
power spectrum that is independent of the matter power 
spectrum. 

Since the shot noise term dominates the polarization 
spectrum, we can understand much of the dependence on 
reionization parameters by approximating the B-mode 
spectrum by its shot noise contribution. For these models 
Cf is well described as a single, wide peak. We will often 
give our results by stating the multipole of this peak, 
^pcak, and its amplitude (Cf ) poak . 

The scale ^ pca k is primarily determined by the bend in 
the averaged one-bubble window function (Wjj), which is 
set by Bcft as described in Section l2~3l More specifically 
^pcak = k c s/DA and since the angular distance varies 
slowly during reionization, the shot noise-dominated B- 
mode power should peak at 

^pcak oc (43) 

with power law behavior as Cf ~ £ 2 on large scales and 
Cf ~ £~ 2 on small scales. Because of the crucial role B ff 
plays in setting the multipole at which the polarization 
spectrum peaks, we present our results in terms of i? c ff 
rather than Rq. 

The amplitude of the peak depends on multiple param- 
eters of the ionization model as we describe in Section |4j 

3.2. Ionization Fraction 

To explore the dependence of the B-modes on x e (z), 
we use a variety of types of ionization histories. 

We assume that x e — for z > Zi and x e = 1 for 
z < Zf. During the epoch of reionization, we require 
< x e < 1. The optical depth places an additional 
constraint on x e (z). The optical depth out to redshift z 
is 

r z (1 + z') 2 
t(z) = cr T n pfi / dz' y ' x e {z'). (44) 
Jq n \ z ) 

We can use the range of observed values for the total 
optical depth r* = t(z > zi) to constrain x e (z) through 
equation (|44[) . Note that the epoch of reionization can 
start no later than sa ~ 12.4 if t* = 0.12, or z% ~ 7.7 if 
r* = 0.06. 

The polarization signal we are concerned with only 
comes from redshifts Zf < z < Zi, so we will de- 
fine the optical depth from the epoch of reionization 
T xe j = — t (zf) , where r(zf) depends only on z/ and 
cosmological parameters that we hold fixed. 

To derive approximate expressions in the following sec- 
tions, we will often assume that the universe is matter 
dominated so that H(z) w B ov / ri m (l + zfl 2 . This is a 
good approximation during the epoch of reionization. 

3.2.1. Constant Ionization Fraction 

The simplest model is a constant value of x e through- 
out the epoch of reionization. For a fixed value of t„, the 
ionization fraction in this model is 



3r r , 



2r 

where tq — <TTn p ,o/( 



(1 + z 4 ) 3/2 - (1 + ^/) 3/2 

^mBo). 



(45) 



7 



Although a model with constant x e is far from being a 
realistic model of reionization, it is a useful simplifying 
assumption for finding the approximate dependence of 
Cf on other parameters of the model. 

3.2.2. Linear Ionization Fraction 

As a more realistic model of reionization we consider an 
ionization fraction that decreases linearly with redshift, 



i — i — I — r 



X e (z) 



1 



z f 



(46) 



For a fixed value of r*, the redshifts Zi and Zf cannot 
be varied independently since they are related through 
equation (|44| . In our linear ionization histories, we fix 
r* and Zf and determine the required value for Zi. The 
upper panel of Figure [5] shows a linear x e (z) for a model 
with t» = 0.12 and zj = 6. 

3.2.3. Maximal B-Mode Ionization History 

If we approximate the H II power spectrum by the shot 
noise term Pfa5 x > an d use the fact that Da(z), Q Tm s(z), 
and t(z) vary slowly with redshift during reionization 
to pull these factors out of the integral in equation ^ . 
then the dependence of the _B-modc amplitude on the 
ionization history is approximately 



C 



dz{l + zf /2 x e {z)[l-x e {z)\. (47) 



We can use variational calculus to determine what form 
of x e (z) maximizes this integral, giving the largest B- 
moclc amplitude in this approximation. 

Fixing the total optical depth r* (and therefore r ro i = 
t* — r(zf)) gives an integral constraint on x e (z), 



/ dz(l + z)V 2 x e (z), 



(48) 



^x e ) + X(l + z)^ 2 x e , (49) 



TO j Zf 

so we can define 

f(z,x e ,x' e ) = (l + zf' 2 x e {l- 

where x' e = dx e /dz and A is a Lagrange multiplier. Since 
this function is independent of x' e , the Euler-Lagrange 
equation is just df/dx e — 0. Solving for x e (z) gives 

x e {z) = 1 + ^(1 + z )-\ (50) 

The ionization fraction given by this equation ap- 
proaches 1/2 at high redshift. For lower values of z, x e 
must stay between by and 1. These limits are satisfied 
for redshifts z > IX] 1 ^ 2 — 1. If z/ satisfies this inequality, 
then x e (z) as given by equation (|50p is between and 1 
during the entire epoch of reionization. If not, then for 
redshifts zf < z < \X\ X / 2 - 1, x e is saturated at either 
or 1 depending on the sign of A. 

The values we assume for t* , Zi, and z f determine the 
sign of A. If A > 0, then x e > 1/2 for all z, so r rc i/ro 
must be greater than it would be with constant x e = 1/2. 
If A < the ionization fraction is always less than 1/2 
and r rc i/ro is less than the constant x e = 1/2 value. We 
can define an optical depth t\ corresponding to the case 
x e (z) = 1/2 



T\{Z % ,Z } ) 



T 



(l + Zi )*/*-(l +Xf ) 



3/2 



(51) 
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Fig. 5. — Upper panel: Ionization history and characteristic 
bubble radius for a linear x e (z) model with Zf = 6, Zi = 18.4, 
t„ = 0.12, and fixed (model 19 in Table Lower panel: 

Ionization histories for three maximal-x e (^) models (models 8-10 
in Table [TJ. The ionization fraction in all three models increases 
to 1 for z < 6. 



so that A in equation (|50|) must have the same sign as 
Trei — T \- For Zf = 6 and z% = 20, t\ ~ 0.1; for Zf = 6 
and Zi = 30, t\ w 0.2. 

After determining the sign of A for particular choices 
of r* , Zi, and z / , the magnitude of A is fixed by the opti- 
cal depth constraint, equation (|48j) . using equation (|50|) 
for x e {z) over the appropriate range of redshifts. The 
resulting ionization histories for three combinations of r* 
and Zi are plotted in the lower panel of Figure [5] 

The maximal-polarization ionization histories are even 
less realistic than those described in the previous sec- 
tions, since they have a significant ionization fraction 
(x e < 1/2) at high redshift, and x e can be small or even 
zero near the end of reionization, if A < 0. These models 
are useful, however, because the B-mode amplitudes as- 
sociated with them provide an upper limit that i?-modes 




x =0.0001 

e 



Q I I I I I I I I I I I I I I I I I I I I I I I L_ 

-8 -6 -4-2 2 
log 10 (n b [Mpc- 3 ]) 

Fig. 6. — Volume averaged bubble radius Ro versus bubble co- 
moving number density for ionization fractions x e = 0.0001, 0.5, 
and 0.9999. Fixed Rq models correspond to horizontal trajecto- 
ries between x e and x e fa 1, and fixed n b models are vertical 
trajectories in this plot. 



from more realistic reionization scenarios must lie below. 

Note that due to the approximations leading to equa- 
tion (|47p . the actual ionization history that maximizes 
the B-mode amplitude will differ slightly from the one 
derived here in equation ([50]) . The results of Monte 
Carlo models that vary x e (z) in each redshift bin sug- 
gest that the true maximal x e (z) is not very different 
from the approximation used here. The B-mode ampli- 
tude of the true maximal-x e model is likely only a few 
percent larger than the model of equation (|50[) , and even 
if the amplitude were much larger it would be only for 
a finely tuned (and unrealistic) ionization history missed 
in the Monte Carlo search. It is reasonable, therefore, to 
take the B polarization produced by models of the ap- 
proximate form derived here as an effective upper limit 
on the amplitude that can be reached by varying x e (z) 
with r* and all other reionization parameters held fixed. 

3.3. Bubble Radius 

We consider two ways of modelling the evolution of 
the bubble size distribution. In one, we assume a fixed 
comoving bubble number density fi& throughout reioniza- 
tion, so that the global ionization fraction increases only 
due to bubble growth. The normalization of the bubble 
radius evolution is set by specifying the comoving char- 
acteristic bubble radius at the end of reionization, R e gj. 
The evolution of R e g in & fixed-fit, model with a linear 
ionization history is shown in the upper panel of FigureO 

The second type of model has a fixed distribution of 
radii, so R c g is independent of z. In this model, the ion- 
ization fraction increases as more ionizing sources col- 
lapse and turn on, while the size of individual bubbles 
remains constant. 

In the radius-number density plane, these models can 
be thought of as trajectories along a horizontal path from 
x e = to x e = 1 for fixed Reg, and vertical trajectories 
for fixed n b (Figure [6]). We assume throughout that a\ n R 



TABLE 1 

MULTIPOLES AND AMPLITUDES OF B POLARIZATION SPECTRA 
PEAKS FOR REIONIZATION MODELS WITH DIFFERENT TYPES OF 
IONIZATION HISTORIES: CONSTANT (C), LINEAR (L), AND FOR 
MAXIMAL S-MODES (M). 



Model 


T* 


Zi 


z f 


fixed 
x e (z) R c s/n b 




RettJ 

[Mpc] 


■^pcak 


(Cg )peak 

[nK 2 ] 


1 


0.12 


20 


6 


C 


both 


1.0 


609 


28 


8836 


2 


0.12 


30 


6 


C 


both 


1.0 


609 


30 


23050 


3 


0.12 


10 


6 


C 


both 


1.0 


609 


32 


45049 


4 


0.06 


9.3 


6 


L 


Rett 


1.0 


609 


24 


422 


5 


0.12 


18.4 


6 


L 


Rett 


1.0 


609 


27 


4325 


6 


0.12 


16.6 


8 


L 


Rett 


1.0 


609 


27 


2944 


7 


0.12 


14.8 


10 


L 


Rett 


1.0 


609 


27 


1600 


8 


0.06 


9.5 


6 


M 


Rett 


1.0 


609 


24 


692 


9 


0.12 


20 


6 


M 


Rett 


1.0 


609 


28 


8989 


10 


0.12 


30 


6 


M 


Rett 


1.0 


609 


31 


27611 


11 


0.12 


10 


6 


M 


Rett 


1.0 


609 


32 


62667 



does not vary with redshift, so the evolution of Rq and 
R is the same as that of i? e ff up to a a\ n ^-dependent 
normalization. 



4. B-MODE PARAMETER DEPENDENCE 

In the previous section, we parameterized the ioniza- 
tion history x e (z) and the effective bubble size R c g(z) 
under several simple schemes that were designed to bring 
out the critical properties of a reionization model for 
the B-mode spectrum. In this section we explore these 
properties and provide scaling relations for the impact of 
reionization parameters on B-modes. 

4.1. Ionization History 

Table [1] gives the .B-mode peak multipoles and ampli- 
tudes for reionization models with the different ionization 
histories described in Section 13.21 For fixed values of t* , 
Zi, Zf, ai n n, and R c g(z), we find that changing the shape 
of the ionization history has only a mild effect on the B- 
mode spectrum compared with the effects of the other 
reionization parameters. 

For example, compare models 1, 5, and 9, which all 
have t* = 0.12, z% ~ 20, Zf — 6, a\ n R = 1.0, and R c g = 
609 comoving Mpc (Rq = 50 Mpc). All of these models 
produce B-modes that peak at £ pca k ~ 28. There is 
about a factor of two difference between the largest peak 
amplitude of these four models (model 9, maximal x e (z) 
with (C^)pcak = 8989 nK 2 ), and the smallest (model 5, 
linear x e {z) with (Cf ) pcak = 4325 nK 2 ). 

The results are similar if we adopt different values for 
the reionization parameters other than x e (z). For ex- 
ample, models 4 and 8, both with t* = 0.06, z.; « 9.5, 
Zf = 6, and the same R c g(z), have the same peak mul- 
tipole, and their B-mode amplitudes differ by less than 
a factor of two between the maximal and linear x e (z) 
models. 

Looking at the entire set of models presented in Ta- 
ble [TJ there is a wide variety of B-mode peak scales 
and amplitudes. However, the largest differences occur 
between models with different choices of r», Zi, Zf, or 
R c b{z). 
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TABLE 2 

MULTIPOLES AND AMPLITUDES OF B POLARIZATION SPECTRA 
PEAKS FOR REIONIZATION MODELS VARYING R B g AND <T lnfl . 



Model 


T* 


Zi 


z f 


fixed 
x e (z) R c s/n b 


°~ln R 


FLre t 
1 t crr , J 

[Mpc] 




\ £. /peak 

[nK 2 ] 


12 


0.12 


30 


6 


C 


both 


1.0 


50 


400 


1928 


13 


0.12 


30 


6 


C 


both 


0.5 


50 


407 


4249 


14 


0.12 


30 


6 


C 


both 


0.2 


50 


450 


6276 


15 


0.12 


30 


6 


c 


both 


1.0 


500 


37 


19057 


16 


0.12 


30 


6 


c 


both 


0.5 


500 


40 


42711 


17 


0.12 


30 


6 


c 


both 


0.2 


500 


45 


63289 


18 


0.06 


9.3 


6 


L 


n b 


1.0 


609 


59 


161 


19 


0.12 


18.4 


6 


L 


n b 


1.0 


609 


72 


1544 


20 


0.12 


16.6 


8 


L 


n b 


1.0 


609 


69 


1092 


21 


0.12 


14.8 


10 


L 


n b 


1.0 


609 


66 


619 



4.2. H II Bubble Size Distribution 

The reionization models in Table [2] and Figure [7] show 
the dependence of £ pca k and (C^P) pca k on H II bubble size. 
Larger values of R e g lead to a smaller multipole for the 
-B-mode peak (see equation (03])) 



^•pcak 



100 



/ 200 Mpc \ 



(52) 



We can see the relation between i? c fj and £ pea k in the 
models in Table [H For example, models 12-14 have dif- 
ferent values of a\ n R but the same i? c ff, and £ pea k is the 
same to within about 10%. As R e s increases, not only 
does ^peak decrease but also the amplitude of the peak 
increases. For a fixed a\ n R it scales as (C^f ) pea k ReS- 

The dependence of the amplitude on the width o~\ n r at 
fixed R e B is much weaker. Note that the B-mode peak 
amplitude for the dominant shot noise contributions can 
be approximated as 



(Cf ) pcak ~ R 3 Q (W 2 (k cS )) 
for fixed i? c ff. Defining 



Kcff 



S R=[- 2 



-7.5a? 



1/4 



«(W 2 R (k cS )), (53) 



(54) 



equation (j3"4")) shows that (W^(fe e ff)) is a function of <7 m r 



only. Evaluating the integral in equation (|34|) . we find 
that e~ 7 - 5a ^ R (W^(k c ff)) is maximized as ai n R — > and 
decreases as ai n R increases, suggesting that (C<f ) peak has 
the same behavior. This happens because the spectrum 
is more sharply peaked as o~\ n r — ► 0, which follows from 
the fact that the window functions (W^(k)) for smaller 
values of o\ n R bend more sharply at k = fc c ff (see Fig- 
ure [3]) ■ This maximum amplitude is estimated to be ap- 
proximately 3 times larger than for the same model with 
cr\nR = 1 (see Figure [7}. 

The models in Table [2] also show this trend with a\ n R, 
and more generally for the range of o~\ n R explored in 
our models, the scaling of the peak amplitude with o\ n R 
at fixed l pea k shows good quantitative agreement with 
equation (|55|) . 

Models 12-17 and the approximate scaling of £ pea k 
and (C^) poa k with R e g assume constant R g throughout 
reionization, but the same general results apply when 
R e g is changing. The only difference is that the constant 
radius must be replaced in the scaling relations by an ef- 
fective radius averaged over the duration of reionization, 
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Fig. 7. — B-mode polarization from inhomogeneous reionization 
with various parameters for the bubble size distribution (models 12- 
17 in Table [2J. The three curves with ^ pcak £S 40 have R eft = 500 
Mpc (solid, thin), and the three with £ pea k ~ 400 have i? e ff = 50 
Mpc (solid, thick), with fixed i? c ff during reionization. For each 
value of R c ff, from the widest spectrum to the narrowest (lowest to 
highest peak amplitude), cr\ n R = 1, 0.5, and 0.2 respectively. For 
all models, the other parameters are Zi = 30, Zf = 6, and r* = 0.12, 
with a constant ionization fraction during reionization. The dashed 
curve shows the predicted B polarization from gravitational lensing 
of i?-modes for our assumed cosmological parameters. 



weighted by the redshift-dependent factors that appear 
in the integral for Cf . Models with fixed-n^ in which 
the radius grows with time are shown in Table [1] (18- 
21). The ^ pea k values for these models indicate that the 
rcdshift-averaged radius for these models is R c s ~ 300 
Mpc (Rq - 25 Mpc). 

For a slowly varying i? c ff and x e , equation J47|) im- 
plies that half the shot noise contributions come from 
0.8 < (1 + z)/(l + Zi) < 1. The _B-mode amplitude 
roughly follows the scalings of the constant R c s model 
with i? e ff evaluated at Z = 0.8(1 + zi) — 1. This expec- 
tation is borne out by models 5 and 19 (see upper panel 
of Figure [5]). More generally for models where R e g de- 
creases with redshift from its final value, as in the fixed- 
fih models, the B-mode peak decreases in amplitude and 
shifts to higher multipoles. 

4.3. Reionization Redshift and Total Optical Depth 

Tableland Figure [8] show models with different values 
of Zi, Zf, and r*. The scalings with these parameters 
can be understood from the shot noise approximation in 
equation (|4T)) 

Cf cx x e (l - x B )[(l + Zl ) 7 / 2 - (1 + z f ) 7 / 2 ] (55) 

if x e is constant during reionization, with x e (T rc j, Zi, zj) 
given by equation (|45|) . Note that the available optical 
depth during the inhomogeneous period r re j sets the ion- 
ization fraction, not the total optical depth r*, and T re i 
is about a factor of 4 larger when = 0.12 than when 
r* = 0.06 if Zf = 6. The ionization fraction between 
these two cases then also differ by roughly a factor of 4, 
but this translates into a somewhat smaller difference in 
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TABLE 3 

MULTIPOLES AND AMPLITUDES OF B POLARIZATION SPECTRA 
PEAKS FOR REIONIZATION MODELS VARYING Zi, Zf, AND T*. 1 5 



Model 


T* 


Zi 


z /' 


fixed 

x e (z) R efl /n b 




[Mpc] 


^-pcak 


irB\ 

\ £ ypcak 

[nK 2 ] 


22 


0.12 


15 


6 


C 


both 


0.5 


93 


194 


962 


23 


0.12 


20 


6 


C 


both 


0.5 


93 


204 


2853 


24 


0.12 


30 


6 


C 


both 


0.5 


93 


217 


7872 


25 


0.12 


10 


6 


C 


both 


0.5 


93 


226 


15154 


26 


0.12 


50 


6 


C 


both 


0.5 


93 


232 


25108 


27 


0.12 


30 


8 


C 


both 


0.5 


93 


218 


6272 


28 


0.12 


30 


10 


C 


both 


0.5 


93 


218 


4025 


29 


0.06 


15 


6 


C 


both 


0.5 


93 


194 


614 


30 


0.06 


50 


6 


C 


both 


0.5 


93 


232 


6792 



the .B-mode amplitudes since Cf depends on j; e (l-i e ) 
and not simply on x e alone. 

The models in Table [3] obey this scaling to good ap- 
proximation. The largest polarization signals are ob- 
tained for models with high Z{ S> Zf where Cf oc (1 + Zj) 2 
(model 26, see also Figure|5]) and external bounds on this 
parameter are the most important to limit B-mode con- 
tributions (see Section [5]). For less extreme cases, the 
scaling of equation (j5l)|) still holds. 

Models 27 and 28 in Table [3] show that ending reion- 
ization earlier leads to smaller B-modes from inhomoge- 
neous reionization. There are two reasons for this: the 
epoch of reionization is shorter, which decreases the term 
in square brackets in equation (|55|) : and increasing Zf 
makes the available optical depth r rc i smaller, which gives 
a smaller value for x e . (Note that this second effect could 
instead act to increase the size of the B polarization if 
x e > 0.5.) One can check that equation (|55|) provides 
a good description of the scaling of B-mode amplitudes 
with Zf to within a few percent accuracy. 

Current observational bounds on the total optical 
depth to reionization allow for a range of values of t* 
that can have a significant effect on the amplitude of the 
B-mode spectrum. Comparing models 26 and 30 in Ta- 
ble [3j we see that increasing r* by a factor of 2 makes the 
amplitude larger by a factor of about 3.7. For the shorter 
epoch of reionization in models 22 and 29, a factor of 2 
increase in r* amplifies the polarization by a factor of 1.6. 
Again, these results are well described by equation (|55"j) 
with x e determined by t». 

In summary, the B-mode amplitude is maximized by 
making the effective bubble radius R e e match the mul- 
tipole of interest, making the duration of the epoch of 
reionization as long as possible, and taking the total opti- 
cal depth to be as large as possible. Narrow distributions 
generate the most signal but the dependence is weak out 

tO CTlni? < 1. 



5. BOUNDS ON REIONIZATION PARAMETERS 

From the previous sections we conclude that the pa- 
rameters that have the greatest effect on the B-mode 
signal from inhomogeneous reionization are R e g and Zi. 
Other parameters can affect the amplitude, but the de- 
gree to which they can do so is limited either on theo- 
retical (x e (z), aini?) or observational grounds (t», Zf). 
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10 100 1000 

1 

Fig. 8. — B-mode polarization spectra from inhomogeneous 
reionization, varying the duration of reionization as in models 22- 
26 in Table l3l (solid curves). From top to bottom, the epoch of 
reionization begins at zi = 50, 40, 30, 20, and 15. For all mod- 
els, the other parameters are Zf = 6, t* = 0.12, o-\ n ji = 0.5, and 
R c ff = 93 Mpc (i?o = 50 Mpc). The models have constant x e and 
R e ff during reionization. The dashed curve shows the predicted 
B-modes from lensing. 




20 25 30 35 40 

Fig. 9. — Contours of Cf at I = 100 in the Zi — R c g plane. 
Other reionization parameters are set to zj = 6, r* = 0.12, and 
^lnii = 0.2, and both x e and R e g are constant during reionization 
in these models. The units of the labeled contours are nK 2 . The 
contour at 2343 nK 2 marks the power in the lensing B-modes at 
£ = 100. The dashed curve in the upper right corner shows the 
comoving radius of a bubble at redshift Zi such that the optical 
depth from that single bubble along a line of sight passing through 
the center of the bubble would be equal to the total optical depth, 
t» = 0.12. The vertical dashed line is at the approximate upper 
limit on zi from the observed shape of the i?-mode reionization 
peak. 
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Figure [9] shows contours of Cf at I = 100 in the Zj — R c s 
plane, with Zf — 6 and r* = 0.12 chosen to maximize 
the amplitude. We take <7\ n R = 0.2 since the amplitude 
at the peak is maximized at small values of <7i n _R. As 
Figure [7] shows, the amplitude on the tails of the spec- 
trum may be higher for larger values of <t\tlR, but the 
difference is no more than a factor of a few. 

The -B-mode amplitude continues to grow as one in- 
creases i? e ff or Zi, so the maximal contributions are de- 
termined by observational or theoretical constraints on 
these quantities. 

5.1. Observational Constraints 

The starting redshift of reionization is limited by the 
shape an d amplitude of t he E-mode peak detected in 
WMAP ()Page et al.ll2006h . The E'-mode peak scale is 
determined by the horizon scale during reionization and 
an early period of reionization shifts the peak to higher 
multipoles than observed. For example, consider a model 
with fixed optical depth t» = 0.09 and constant ioniza- 
tion fraction out to Zi = 30. In this model, the suppres- 
sion of power at £ = 3, where the weight of the detec- 
tion lies, is a factor of 2 larger than in an instantaneous 
reionization model with r» = 0.06, while the enhance- 
ment of power at £ — 10 is a factor of 4 larger than in a 
model with t* = 0-12. Thi s is consistent with the find- 
ings of iSpergel et al.l (J2006) for constant x e models out 
to Zi = 25 and implies that Zi < 30 at approximately the 
2a level. We show this bound as the vertical dashed line 
in Figure [9] 

With this bound, models that have B-mode power 
in excess of the lensing signal at £ — 100 have i? e ff ^ 
70 Mpc. Bubbles with large radii at high redshift also 
contain a large optical depth. The optical depth from a 
fully ionized bubble at redshift z with radius R along a 
line of sight that passes through the center of the bubble 
is 

n = o- T n p , (l + z) 2 2R. (56) 

This single bubble optical depth can exceed the total 
r* if such bubbles have sufficiently low number densities. 
Despite being rare, they still produce an observable effect 
on the temperature power spectrum. In the direction 
of such a bubble, the temperature power spectrum will 
be suppressed by an additional factor of e~ 2Tb . These 
bubbles would produce a power spectrum of the acoustic 
peaks which varied from field to field and experiment 
to experiment. No such effect has been seen and the 
differences between measurements can be attributed to 
the ~ 5-10% calibration errors. 

Such an effect is also limited by the fiel d-to-field vari- 
ance of the WMAP power spectrum. IHansen et al.l 
(2004) found no evidence for excess variance beyond the 
10% sample variance of the amplitude of the first peak 
when measured in different 19°-radius patches on the sky. 

As a conservative limit, we require that the optical 
depth from a single bubble should be smaller than the 
total optical depth to reionization r* = 0.12 or a 20% 
variation in the power in the acoustic peaks on and off 
of a bubble. For the radius, this constraint becomes 

fl < 14 0Mpo(^)(l±£)" 2 . (57, 

Since the typical size of a bubble that contributes to the 
B polarization from inhomogeneous reionization is R c s, 



we apply this constraint to the effective radius. This 
bound is shown for Tb — 0.12 as the dashed curve in 
Figure O Combined with equation (f52"| . this also implies 
a bound on the peak multipole of 




Therefore even if the bound on Zi is relaxed the B-mode 
contributions at £ ~ 100 would be in the steeply falling 
white noise tail of the spectrum. Correspondingly the 
maximal contributions for models in Figure[5]that satisfy 
this constraint occur for z% < 25. 

This constraint on R c g scales linearly wi th the optical 
depth Tb, so future work similar to that of IHansen et alJ 
(|2004h could potentially strengthen this bound. For ex- 
ample, a bound of Tb < 0.06 would exclude models with 
B-modes above the lensing signal at Zi = 30. 

These bounds can potentially be evaded if the width 
of the distribution o\ n u is so wide that bubbles with 
R ~ B G ff represent such a small fraction of the total that 
a typical line of sight would never intersect such a bubble. 
However, if such bubbles are the source of the B-modcs, 
then there would be a large fraction of sky that would 
be free of this signal. These clean regions of sky could 
then be used to measure gravitational wave B-modes. To 
quantify these considerations, we calculate the covering 
fraction of the sky in bubbles with R ~ R e g as a function 

Of CTlnfl- 

We define the covering fraction as follows. Consider a 
radius R p such that p% of the power in the one-bubble 
shot noise term comes from bubbles with radii R > R p . 
The median bubble radius that contributes to B-modes 
is B50 , which is within 5-20% of R c s for the range of a\ n r 
that we consider. We select bubbles with R ~ B e ff that 
contribute half the total power by requiring bubble radii 
to satisfy R 75 < R < B25. For a\ nR = 0.2, R 75 = 0.9B cff 
and R25 = 1.2B c ff; for cr\ n R = 1, B75 = 0.8B c ff and 
B25 = 1.6B e ff. 

We also choose a representative range in redshift for the 
covering fraction. Since most of the polarization signal 
comes from high-redshift bubbles, we select redshifts zi — 
Az < z < zi at which bubbles contribute half the total 
B-mode power. In the reionization models we consider, 
Az varies between about 0.15(1 + ^) and 0.25(1 + z,-), so 
we approximate the width in redshift as Az w 0.2(1 + Zi) 
(see Section 14^2^1 . Since bubbles with radii between B75 
and B25 contribute half the total power at any given 
redshift, combining the restrictions on bubble radius and 
redshift leaves a sample of bubbles that produce a quarter 
of the B-mode power from inhomogeneous reionization. 

In a spherical shell at redshift Zi with width Az, the 
number of bubbles per solid angle with radius R is 

and the solid angle subtended by a bubble of radius R at 
distance Da is ttR 2 /D a . Then the covering fraction of 
bubbles that meet the criteria for radius and redshift is 

f R25 , dN irR 2 , . 

fcovci = Jn 75 dR dRdn^\ m 

= 3(1 + Zj)x e c 3al„/2 

40H( Zl )R cS 
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Fig. 10. — The covering fraction of bubbles with radii R ~ R e g 
at redshifts z < Zi that are responsible for about 1/4 of the total 
shot noise power. The covering fraction depends on i? e fj, Zi, and 
indicated in the axis label. 



erf 



V V2ai n 



R 



erf 



2cr lni?. + lnr 75 



\f2a\n 



R 



where we define r p = R p /R c g and assume Az = 0.2(1 + 
Zi) and i e C 1. Note that this is likely to be an overesti- 
mate as the covering fraction approaches unity since we 
assume that the bubbles are not overlapping along lines 
of sight. 

The covering fraction for z% — 30, R c g = 50 Mpc, and 
x e = 0.1 is shown in Figure flOl It becomes negligible 
for distributions with a\ n R > 1. The covering fraction 
decreases for larger z,, larger i? e ff, or smaller x e . 

5.2. Theoretical Expectations 

Observational constraints still allow models with large 
B-mode contributions, approaching and in some cases 
exceeding the gravitational lens contributions (see Fig- 
ure [5]). However they require R e g > 100 Mpc at z% ~ 20- 
30. These models are currently disfavored theoretically 
as plausible reionization scenarios. 

Normal population II star formation is unlikely to be 
suffici ent to cause signific ant ionization above z ~ 15 
(see lBarkana fe Loebll2001l . and references therein). One 
possibility for obtaining substantial ionization at high 
redshift is a "double reionization" model where a change 
in the ionizing sources to say metal free population III 
stars ca uses non-monotonic behavior in the ionization 
fraction dCenl 120031 : IWvithe fc Loebl[200l . Such models 
were studied extensively following the indications from 
the first ye ar WMAP data of a high t otal optical depth. 

Recently, iFurlanetto fe Loebl ((2005) have argued that 
such models are unlikely to cause a non-monotonic be- 
havior. This would rule out models that have signifi- 
cant ionization at high redshifts compensated by a low 
value of x e at intermediate redshifts to keep the total 
optical depth within the observed limits, such as the 
maxi mal .B-mode ionizat i on his tories we consider. How- 
ever, IFurlanetto fe Loebl (|2005l ) find that models with 



extended periods of ionization in which x e is significant 
over a wide range of redshifts out to Zj < 20 are not 
difficult to find. Recall that such models would produce 
B-modes that are within a factor of a few of maximal 
at a given Z{. Ho wever given the low as currently indi- 
cated by WMAP (|Spergel et all 120061) . extended^ reion- 
ization at high redshift beco mes even more difficult to 
achieve (|Alvarez et al.l 12006) . Recent simulations have 
nonetheless found that substantial ionization at z = 15- 
20 is still possible with small-mass sources (|Iliev et al.l 
l2006h . 

Even for Zi < 20 substantial .B-modes are possible but 
only if the effective bubble radius R e s > 100 Mpc. The 
clustering of sources could create H II regions with co- 
moving radii of 10-30 Mpc near the end of reionization, 
larger than might be expected for a bub ble around a sin- 
gle source galaxy (jWvithe fe Loebll2005D . Note also that 
i? e ff is much greater than the typical radius R of a bub- 
ble if the distribution is wide. However, dense absorbers 
such as Lyman limit systems that are found to be as- 
sociated with ga laxies in simulations could restrict the 
size o f bubbles (|Gnedin fe Fan! [200l iKohler fc Giiediil 
2007). Furthermore, near the beginning of reionization 
the existing bubbles should be much smaller in radius. 

To make these c onsiderations concrete, w e examine the 
analytic model of IFurlanetto et all (|2004f ) which makes 
specific predictions for the bub ble distribu t ion w hen 
coupled with the simulations of IZahn et"aT1 (|2007l ) for 
the ionization history. Here reionization occurs fairly 
promptly and at low redshift. In this model, the peak 
radius of the size distribution is of order R ~ 1 Mpc (co- 
moving) early in reionization and increases to R ~ 10-10 
Mpc by the end of reionization (jFurlanetto et alj 120061) . 
The shape of the distribution shows some skewness to- 
ward smaller bubbles but can be crudely modelled by a 
log normal with a width that increases with redshift from 
fini?, ~ 0.5 to crinfl ~ 1. For simplicity we approximate 
this as a constant a\ n R = 0.7. Both the ionization history 
and the bubble distribution in this model act to lower the 
B-mode contributions as shown in Figure QTJ They are 
several orders of magnitude below the lensing .B-modes 
and would be relevant only if the polarization field were 
cleaned of lensing on the whole sky. This model can also 
be considered near the minimal expected contribution in 
that the total optical depth is quite low (r* = 0.055) and 
the duration of reionization quite short. 

The model in Figure [TT] illustrates some exceptions 
to our earlier conclusions. Because x e and i? c ff change 
rapidly with redshift, the redshift-averaged effective ra- 
dius is not necessarily well approximated by R c g(0.8zi) 
as described in Section 14.21 The location of the peak 
at ^p C ak ~ 400 indicates that the main contribution is 
from bubbles with radii R ~ 50 Mpc, and R e g does not 
reach this size until around z = 7. Related to this is the 
fact that in this model, nearly all of the B-mode power 
comes from bubbles at z < 7, at the very end of reion- 
ization. In the models we examined where reionization 
is a more gradual process, a large fraction of the power 
comes instead from the highest-redshift bubbles. 

Our model assumes that high-density regions that 
host ionizing sources are ionized before areas with lower 
density, co nsistent wi t h recent simulations of reioniza- 
tion ( e.g.. ICenl 120031 : ISokasian et all 12004 llliev et all 
120061 IGnedin fc Fanfl2006t IZahn et alJ 12007^ Mod- 
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Fig. 11. — Ionization history and bubble size (upper panel) and 
B polarization from inhomogeneous reionization (lower panel) for 
a model inten ded to appr o ximat e the results of the reionization 
simulations of IZahn et al.l :|2IM1",':. The width of the bubble size 
distribution is chosen to be 0"i n ij = 0.7. As in Figure [4] the tens- 
ing S-modes (dashed) and gravitational wave S-modes (dotted) 
are shown for comparison. The components of the inhomogeneous 
reionization power spectrum are also plotted as in Figure [4] The 
total optical depth for this model is r* = 0.055 and the duration 
of reionization is fairly short so this model can be taken as nearly 
the minimal B-modes from inhomogeneous reionization. 



els in which reionization proceeds in the reverse order 
have also been proposed, and such a scenario could 
produ ce different .B-mode spectr a from those studied 
here (Miralda-Esc ude et alj [2000h . If reionization does 
begin in voids, we would expect the two-bubble contri- 
bution to the B polarization to be affected the most. As 
long as the H II regions in such a scenario could still be 
characterized as roughly spherical ionized regions with a 
size distribution that is approximately log normal, our 
model should still be sufficient to describe the shot noise 
contribution to the B-modes. 



6. DISCUSSION 

We have presented a comprehensive study of the gen- 
eration of 5-mode polarization by inhomogeneous reion- 
ization with an emphasis on the phenomenological prop- 
erties required to generate a substantial contamination 
for degree scale gravitational wave studies. We base our 
study on a general parameterization of inhomogeneous 
reionization in terms of three functions of redshift: the 
mean ionization fraction, the effective bubble radius, and 
the width of a log normal radius distribution. 

We find that these -B-modes are only important when 
the ionized bubbles are sufficiently rare that the power 
spectrum is dominated by the shot noise of the bubbles. 
In such models, the -B-modes are maximized by taking 
an effective bubble radius as large as possible, up to a 
few hundred Mpc, at as high a redshift as possible for 
a total duration of reionization that is as long as possi- 
ble, while still remaining consistent with constraints on 
the total optical depth. The width of the bubble distri- 
bution is less critical but a wider distribution can make 
the effective bubble radius much larger than the typical 
bubble. The details of the ionization history are also less 
relevant so long as a finite ionization fraction remains at 
high redshift. We provide useful scaling relations for es- 
timating the 5-mode contributions for a wide variety of 
scenarios. 

These conditions are constrained both observationally 
and by theoretical expectations. A substantial ionization 
fraction at high redshift would violate the shape of the 
observed -B-mode power spectrum which sets a maximum 
redshift of z ~ 30. A large bubble radius would produce 
an large optical depth through its center and hence in- 
homogeneous obscuration of the acoustic peaks. Com- 
bined these conditions still allow power at I = 100 that 
is ~ 5 times the power produced by gravitational lensing 
or ~ 0.1 fiK in the polarization field. These amplitudes 
can be reached for example by ~ 140 Mpc bubbles at 
z ~ 30. This level is also comparable to the maximum 
allowed -B-modes from gravitational waves at £ = 100. 
On the other hand, not even these models would pre- 
dict -B-modes at I < 10 that would mask the maximal 
gravitational wave signal from reionization. 

Current theoretical models and simulations would not 
predict the existence of such large ionization bubbles at 
such high redshifts. In these models bubbles grow to a 
few tens of Mpc only toward the end of a fairly prompt 
reionization. Nonetheless, even for conditions where the 
.B-modes do not exceed the lensing signal, they can still 
be the leading order cosmological contaminant to the 
gravitational wave signal. 

Even without de-lensing of the polarization, the well- 
defined power spectrum from lensing allows its contri- 
butions to be statistically subtracted, ultimately to the 
few percent level with a full-sky cosmic variance limited 
experiment. The uncertain form of the inhomogeneous 
reionization .B-modes does not permit statistical subtrac- 
tion. With r* = 0.12, ~ 20 Mpc bubbles at z ~ 15 near 
the end of reionization would produce -B-modes that are 
1% of the power in lensing at £ = 100 or 10% of the 
polarization field. Likewise with a lower total optical 
depth of t* = 0.06, ~ 40 Mpc bubbles at z ~ 10 would 
produce the same effect. Such contributions would also 
be relevant for the percent precision measurements of 
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the lensing _S-modes on smaller scales expected from the 
next generation of ground based experiments. 
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